Embryonic and juvenile snakes (Natrix maura, Linnaeus 1758) compensate for high elevation hypoxia via shifts in cardiovascular physiology and metabolism.

The colonization of novel environments requires a favorable response to conditions never, or rarely, encountered in recent evolutionary history. For example, populations colonizing upslope habitats must cope with lower atmospheric pressure at elevation, and thus reduced oxygen availability. The embryo stage in oviparous organisms is particularly susceptible, given its lack of mobility and limited gas exchange via diffusion through the eggshell and membranes. Especially little is known about responses of Lepidosaurian reptiles to reduced oxygen availability. To test the role of physiological plasticity during early development in response to high elevation hypoxia, we performed a transplant experiment with the viperine snake (Natrix maura, Linnaeus 1758). We maintained gravid females originating from low elevation populations (432 m above sea level [ASL]-normoxia) at both the elevation of origin and high elevation (2877 m ASL-extreme high elevation hypoxia; approximately 72% oxygen availability relative to sea level), then incubated egg clutches at both low and high elevation. Regardless of maternal exposure to hypoxia during gestation, embryos incubated at extreme high elevation exhibited altered developmental trajectories of cardiovascular function and metabolism across the incubation period, including a reduction in late-development egg mass. This physiological response may have contributed to the maintenance of similar incubation duration, hatching success, and hatchling body size compared to embryos incubated at low elevation. Nevertheless, after being maintained in hypoxia, juveniles exhibit reduced carbon dioxide production relative to oxygen consumption, suggesting altered energy pathways compared to juveniles maintained in normoxia. These findings highlight the role of physiological plasticity in maintaining rates of survival and fitness-relevant phenotypes in novel environments.

Lizards from warm and declining populations are born with extremely short telomeres.

Aging is the price to pay for acquiring and processing energy through cellular activity and life history productivity. Climate warming can exacerbate the inherent pace of aging, as illustrated by a faster erosion of protective telomere DNA sequences. This biomarker integrates individual pace of life and parental effects through the germline, but whether intra- and intergenerational telomere dynamics underlies population trends remains an open question. Here, we investigated the covariation between life history, telomere length (TL), and extinction risk among three age classes in a cold-adapted ectotherm (Zootoca vivipara) facing warming-induced extirpations in its distribution limits. TL followed the same threshold relationships with population extinction risk at birth, maturity, and adulthood, suggesting intergenerational accumulation of accelerated aging rate in declining populations. In dwindling populations, most neonates inherited already short telomeres, suggesting they were born physiologically old and unlikely to reach recruitment. At adulthood, TL further explained females' reproductive performance, switching from an index of individual quality in stable populations to a biomarker of reproductive costs in those close to extirpation. We compiled these results to propose the aging loop hypothesis and conceptualize how climate-driven telomere shortening in ectotherms may accumulate across generations and generate tipping points before local extirpation.

Viability, behavior, and color expression in the offspring of matings between common wall lizard Podarcis muralis color morphs.

Color polymorphisms are widely studied to identify the mechanisms responsible for the origin and maintenance of phenotypic variability in nature. Two of the mechanisms of balancing selection currently thought to explain the long-term persistence of polymorphisms are the evolution of alternative phenotypic optima through correlational selection on suites of traits including color and heterosis. Both of these mechanisms can generate differences in offspring viability and fitness arising from different morph combinations. Here, we examined the effect of parental morph combination on fertilization success, embryonic viability, newborn quality, antipredator, and foraging behavior, as well as inter-annual survival by conducting controlled matings in a polymorphic lacertid Podarcis muralis, where color morphs are frequently assumed to reflect alternative phenotypic optima (e.g., alternative reproductive strategies). Juveniles were kept in outdoor tubs for a year in order to study inter-annual growth, survival, and morph inheritance. In agreement with a previous genome-wide association analysis, morph frequencies in the year-old juveniles matched the frequencies expected if orange and yellow expressions depended on recessive homozygosity at 2 separate loci. Our findings also agree with previous literature reporting higher reproductive output of heavy females and the higher overall viability of heavy newborn lizards, but we found no evidence for the existence of alternative breeding investment strategies in female morphs, or morph-combination effects on offspring viability and behavior. We conclude that inter-morph breeding remains entirely viable and genetic incompatibilities are of little significance for the maintenance of discrete color morphs in P. muralis from the Pyrenees.

Impacts of temperature on O2 consumption of the Pyrenean brook newt (Calotriton asper) from populations along an elevational gradient.

Global warming impacts biodiversity worldwide, leading to species' adaptation, migration, or extinction. The population's persistence depends on the maintenance of essential activities, which is notably driven by phenotypic adaptation to local environments. Metabolic rate - that increases with temperature in ectotherms - is a key physiological proxy for the energy available to fuel individuals' activities. Cold-adapted ectotherms can exhibit a higher resting metabolism than warm-adapted ones to maintain functionality at higher elevations or latitudes, known as the metabolic cold-adaptation hypothesis. How climate change will affect metabolism in species inhabiting contrasting climates (cold or warm) is still a debate. Therefore, it is of high interest to assess the pace of metabolic responses to global warming among populations adapted to highly different baseline climatic conditions. Here, we conducted a physiological experiment in the endemic Pyrenean brook newt (Calotriton asper). We measured a proxy of standard metabolic rate (SMR) along a temperature gradient in individuals sampled among 6 populations located from 550 to 2189 m a.s.l. We demonstrated that SMR increased with temperature, but significantly diverged depending on populations' origins. The baseline and the slope of the relationship between SMR and temperature were both higher for high-elevation populations than for low-elevation populations. We discussed the stronger metabolic response observed in high-elevation populations suggesting a drop of performance in essential life activities for these individuals under current climate change. With the increase of metabolism as the climate warms, the metabolic-cold adaptation strategy selected in the past could compromise the sustainability of cold-adapted populations if short-term evolutionary responses do not allow to offset this evolutionary legacy.

Metal(loid) pollution, not urbanisation nor parasites predicts low body condition in a wetland bioindicator snake.

Urban ecosystems and remnant habitat 'islands' therein, provide important strongholds for many wildlife species including those of conservation significance. However, the persistence of these habitats can be undermined if their structure and function are too severely disrupted. Urban wetlands, specifically, are usually degraded by a monoculture of invasive vegetation, disrupted hydrology, and chronic-contamination from a suite of anthropogenic pollutants. Top predators-as bioindicators-can be used to assess and monitor the health of these ecosystems. We measured eight health parameters (e.g., parasites, wounds and scars, tail loss and body condition) in a wetland top predator, the western tiger snake, Notechis scutatus occidentalis. For three years, snakes were sampled across four wetlands along an urban gradient. For each site, we used GIS software to measure the area of different landscapes and calculate an urbanisation-landscape score. Previously published research on snake contamination informed our calculations of a metal-pollution index for each site. We used generalised linear mixed models to assess the relationship between all health parameters and site variables. We found the metal-pollution index to have the most significant association with poor body condition. Although parasitism, tail loss and wounds differed among sites, none of these parameters influenced body condition. Additionally, the suite of health parameters suggested differing health status among sites; however, our measure of contemporary landscape urbanisation was never a significant predictor variable. Our results suggest that the health of wetland predators surrounding a rapidly growing city may be offset by higher levels of environmental pollution.

Population Genomics of Wall Lizards Reflects the Dynamic History of the Mediterranean Basin.

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an "Out-of-Italy" expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

Thermal-metabolic phenotypes of the lizard Podarcis muralis differ across elevation, but converge in high-elevation hypoxia.

In response to a warming climate, many montane species are shifting upslope to track the emergence of preferred temperatures. Characterizing patterns of variation in metabolic, physiological and thermal traits along an elevational gradient, and the plastic potential of these traits, is necessary to understand current and future responses to abiotic constraints at high elevations, including limited oxygen availability. We performed a transplant experiment with the upslope-colonizing common wall lizard (Podarcis muralis) in which we measured nine aspects of thermal physiology and aerobic capacity in lizards from replicate low- (400 m above sea level, ASL) and high-elevation (1700 m ASL) populations. We first measured traits at their elevation of origin and then transplanted half of each group to extreme high elevation (2900 m ASL; above the current elevational range limit of this species), where oxygen availability is reduced by ∼25% relative to sea level. After 3 weeks of acclimation, we again measured these traits in both the transplanted and control groups. The multivariate thermal-metabolic phenotypes of lizards originating from different elevations differed clearly when measured at the elevation of origin. For example, high-elevation lizards are more heat tolerant than their low-elevation counterparts (counter-gradient variation). Yet, these phenotypes converged after exposure to reduced oxygen availability at extreme high elevation, suggesting limited plastic responses under this novel constraint. Our results suggest that high-elevation populations are well suited to their oxygen environments, but that plasticity in the thermal-metabolic phenotype does not pre-adapt these populations to colonize more hypoxic environments at higher elevations.

Bioindicator snake shows genomic signatures of natural and anthropogenic barriers to gene flow.

Urbanisation alters landscapes, introduces wildlife to novel stressors, and fragments habitats into remnant 'islands'. Within these islands, isolated wildlife populations can experience genetic drift and subsequently suffer from inbreeding depression and reduced adaptive potential. The Western tiger snake (Notechis scutatus occidentalis) is a predator of wetlands in the Swan Coastal Plain, a unique bioregion that has suffered substantial degradation through the development of the city of Perth, Western Australia. Within the urban matrix, tiger snakes now only persist in a handful of wetlands where they are known to bioaccumulate a suite of contaminants, and have recently been suggested as a relevant bioindicator of ecosystem health. Here, we used genome-wide single nucleotide polymorphism (SNP) data to explore the contemporary population genomics of seven tiger snake populations across the urban matrix. Specifically, we used population genomic structure and diversity, effective population sizes (Ne), and heterozygosity-fitness correlations to assess fitness of each population with respect to urbanisation. We found that population genomic structure was strongest across the northern and southern sides of a major river system, with the northern cluster of populations exhibiting lower heterozygosities than the southern cluster, likely due to a lack of historical gene flow. We also observed an increasing signal of inbreeding and genetic drift with increasing geographic isolation due to urbanisation. Effective population sizes (Ne) at most sites were small (< 100), with Ne appearing to reflect the area of available habitat rather than the degree of adjacent urbanisation. This suggests that ecosystem management and restoration may be the best method to buffer the further loss of genetic diversity in urban wetlands. If tiger snake populations continue to decline in urban areas, our results provide a baseline measure of genomic diversity, as well as highlighting which 'islands' of habitat are most in need of management and protection.

PLASMA BIOCHEMISTRY PROFILES OF WILD WESTERN TIGER SNAKES (NOTECHIS SCUTATUS OCCIDENTALIS) BEFORE AND AFTER SIX MONTHS OF CAPTIVITY.

Urban wildlife often suffer poorer health than their counterparts living in more pristine environments due to exposure to anthropogenic stressors such as habitat degradation and environmental contamination. As a result, the health of urban versus nonurban snakes might be assessed by differences in their plasma biochemistries. We compared the plasma profiles of western tiger snakes (Notechis scutatus occidentalis) from a heavily urbanized wetland and a natural, nonurbanized wetland. Despite the urbanized snakes having lower body mass index, we found no significant difference between the plasma profiles of the two populations. We collected snakes from each population and kept them in captivity for 6 mo, providing them with stable conditions, uncontaminated (exempt from heavy metals and pesticides) food and water, and lowered parasite intensity in an attempt to promote better health through depuration. After captivity, snakes experienced a significant improvement in body mass index and significant changes in their plasma profiles. Snakes from the natural wetland initially had more variation of DNA damage; mean concentration of DNA damage in all snakes slightly decreased, but not significantly, after captivity. We present the plasma biochemistry profiles from western tiger snakes both before and after captivity and suggest a period of removal from natural stressors via captivity may offer a more reliable result of how plasma profiles of healthy animals might appear.

No evidence for differential sociosexual behavior and space use in the color morphs of the European common wall lizard (Podarcis muralis).

Explaining the evolutionary origin and maintenance of color polymorphisms is a major challenge in evolutionary biology. Such polymorphisms are commonly thought to reflect the existence of alternative behavioral or life-history strategies under negative frequency-dependent selection. The European common wall lizard Podarcis muralis exhibits a striking ventral color polymorphism that has been intensely studied and is often assumed to reflect alternative reproductive strategies, similar to the iconic "rock-paper-scissors" system described in the North American lizard Uta stansburiana. However, available studies so far have ignored central aspects in the behavioral ecology of this species that are crucial to assess the existence of alternative reproductive strategies. Here, we try to fill this gap by studying the social behavior, space use, and reproductive performance of lizards showing different color morphs, both in a free-ranging population from the eastern Pyrenees and in ten experimental mesocosm enclosures. In the natural population, we found no differences between morphs in site fidelity, space use, or male-female spatial overlap. Likewise, color morph was irrelevant to sociosexual behavior, space use, and reproductive success within experimental enclosures. Our results contradict the commonly held hypothesis that P. muralis morphs reflect alternative behavioral strategies, and suggest that we should instead turn our attention to alternative functional explanations.

Comment on "Forest microclimate dynamics drive plant responses to warming".

Zellweger et al (Reports, 15 May 2020, p. 772) claimed that forest plant communities' response to global warming is primarily controlled by microclimate dynamics. We show that community thermophilization is poorly explained by the underlying components of microclimate, and that global warming primarily controls the climatic lag of plant communities. Deconstructing the underlying components of microclimate provides insights for managers.

High-elevation hypoxia impacts perinatal physiology and performance in a potential montane colonizer.

Climate change is generating range shifts in many organisms, notably along the elevational gradient in mountainous environments. However, moving up in elevation exposes organisms to lower oxygen availability, which may reduce the successful reproduction and development of oviparous organisms. To test this possibility in an upward-colonizing species, we artificially incubated developing embryos of the viperine snake (Natrix maura) using a split-clutch design, in conditions of extreme high elevation (hypoxia at 2877 m above sea level; 72% sea-level equivalent O2 availability) or low elevation (control group; i.e. normoxia at 436 m above sea level). Hatching success did not differ between the two treatments. Embryos developing at extreme high elevation had higher heart rates and hatched earlier, resulting in hatchlings that were smaller in body size and slower swimmers compared to their siblings incubated at lower elevation. Furthermore, post-hatching reciprocal transplant of juveniles showed that snakes which developed at extreme high elevation, when transferred back to low elevation, did not recover full performance compared to their siblings from the low elevation incubation treatment. These results suggest that incubation at extreme high elevation, including the effects of hypoxia, will not prevent oviparous ectotherms from producing viable young, but may pose significant physiological challenges on developing offspring in ovo. These early-life performance limitations imposed by extreme high elevation could have negative consequences on adult phenotypes, including on fitness-related traits.

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees.

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present-day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present-day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex-biased dispersal. This work highlights the importance of integrating past evolutionary processes and present-day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.

Vertical Transmission of a Nematode from Female Lizards to the Brains of Their Offspring.

Parasites have evolved a diversity of lifestyles that exploit the biology of their hosts. Some nematodes that parasitize mammals pass via the placenta or milk from one host to another. Similar cases of vertical transmission have never been reported in avian and nonavian reptiles, suggesting that egg laying may constrain the means of parasite transmission. However, here we report the first incidence of transovarial transmission of a previously undescribed nematode in an egg-laying amniote, the common wall lizard (Podarcis muralis). Nematodes enter the developing brain from the female ovary early in embryonic development. Infected lizard embryos develop normally and hatch with nematodes residing in their braincase. We present a morphological and molecular phylogenetic characterization of the nematode and suggest that particular features of lizard biology that are absent from birds and turtles facilitated the evolutionary origin of this novel life history.

A simple ATAC-seq protocol for population epigenetics.

We describe here a protocol for the generation of sequence-ready libraries for population epigenomics studies. The protocol is a streamlined version of the Assay for transposase accessible chromatin with high-throughput sequencing (ATAC-seq) that provides a positive display of accessible, presumably euchromatic regions. The protocol is straightforward and can be used with small individuals such as daphnia and schistosome worms, and probably many other biological samples of comparable size, and it requires little molecular biology handling expertise.

Investigating the role of urbanisation, wetlands and climatic conditions in nematode parasitism in a large Australian elapid snake.

Tiger snakes (Notechis scutatus) in wetlands of South-West Western Australia (SW WA) are commonly parasitised by the nematode Ophidascaris pyrrhus. Host-parasite interactions are complex and can potentially be impacted by factors such as urbanisation or climate. We assessed whether urbanisation, distance to wetland sites, and climatic factors have influenced parasitism in tiger snakes from specimens collected over the last century. We dissected 91 museum specimens of tiger snakes across SW WA and counted gastrointestinal nematodes. Binomial generalised linear modelling, with presence/absence of nematodes as a response variable, was used to determine which factors were driving infection. Model selection using AICc values showed that proximity to wetlands, rainfall and topographic wetness were most strongly associated with the probability of infection of snakes by nematodes. We also found a slight positive correlation between nematode abundance and annual mean maximum temperature. We found no significant influence of distance to urban centre on nematode burdens; however, our results suggest that water-related variables are a key driver of nematode parasitism in tiger snakes in SW WA. We also suggest that urbanisation is still of interest as its role in wetland and climate modification may increase parasitism in wetland snakes.

Isolation and characterization of fourteen polymorphic microsatellite markers in the viperine snake Natrix maura.

Nineteen polymorphic microsatellite loci were identified and developed for Natrix maura. Polymorphism was assessed for 120 individuals sampled across four sampling sites from the French Pyrenees Mountains. The number of alleles per locus ranged from 3 to 15, and expected heterozygosity per locus ranged from 0.227 to 0.863. We tested for deviation from Hardy-Weinberg equilibrium and linkage disequilibrium and assessed the presence of null alleles for all loci, resulting in a selection of 14 high-quality polymorphic markers. These markers will be extremely useful in identifying fine-scale genetic structures and providing insight into conservation management plans of this species.

Transplanting gravid lizards to high elevation alters maternal and embryonic oxygen physiology, but not reproductive success or hatchling phenotype.

Increased global temperatures have opened previously inhospitable habitats, such as at higher elevations. However, the reduction of oxygen partial pressure with increase in elevation represents an important physiological constraint that may limit colonization of such habitats, even if the thermal niche is appropriate. To test the mechanisms underlying the response to ecologically relevant levels of hypoxia, we performed a translocation experiment with the common wall lizard (Podarcis muralis), a widespread European lizard amenable to establishing populations outside its natural range. We investigated the impacts of hypoxia on the oxygen physiology and reproductive output of gravid common wall lizards and the subsequent development and morphology of their offspring. Lowland females transplanted to high elevations increased their haematocrit and haemoglobin concentration within days and maintained routine metabolism compared with lizards kept at native elevations. However, transplanted lizards suffered from increased reactive oxygen metabolite production near the oviposition date, suggesting a cost of reproduction at high elevation. Transplanted females and females native to different elevations did not differ in reproductive output (clutch size, egg mass, relative clutch mass or embryonic stage at oviposition) or in post-oviposition body condition. Developing embryos reduced heart rates and prolonged incubation times at high elevations within the native range and at extreme high elevations beyond the current range, but this reduced oxygen availability did not affect metabolic rate, hatching success or hatchling size. These results suggest that this opportunistic colonizer is capable of successfully responding to novel environmental constraints in these important life-history stages.

Pattern matters: Snakes exhibiting triangular and diamond-shaped skin patterns modulate electrophysiological activity in human visual cortex.

The neural and perceptual mechanisms that support the efficient visual detection of snakes in humans are still not fully understood. According to the Snake Detection Theory, selection pressures posed by snakes on early primates have shaped the development of the visual system. Previous studies in humans have investigated early visual electrophysiological activity in response to snake images vs. various alternative dangerous or non-dangerous stimuli. These studies have shown that the Early Posterior Negativity (EPN) component is selectively elicited by snake or snake-like images. Recent findings yielded the complementary/alternative hypothesis that early humans (and possibly other primates) evolved an aversion especially for potentially harmful triangular shapes, such as teeth, claws or spikes. In the present study we investigated the effect of triangular and diamond-shaped patterns in snake skins on the ERP correlates of visual processing in humans. In the first experiment, we employed pictures of snakes displaying either triangular/diamond-shaped patterns or no particular pattern on their skins, and pictures of frogs as control. Participants observed a random visual presentation of these pictures. Consistent with previous studies, snakes elicited an enhanced negativity between 225 and 300 ms (EPN) compared to frogs. However, snakes featuring triangular/diamond-shaped patterns on their skin produced an enhanced EPN compared to the snakes that did not display such patterns. In a second experiment we used pictures displaying only skin patterns of snakes and frogs. Results from the second experiment confirmed the results of the first experiment, suggesting that triangular snake-skin patterns modulate the activity in human visual cortex. Taken together, our results constitute an important contribution to the snake detection theory.

Lizards at the Peak: Physiological Plasticity Does Not Maintain Performance in Lizards Transplanted to High Altitude.

Warming climates are facilitating the range expansion of many taxa to habitats that were formerly thermally inhospitable, including to higher latitudes and elevations. The potential for such colonization, however, varies widely among taxa. Because environmental factors may interact to affect colonization potential, an understanding of underlying physiological and behavioral mechanisms is necessary to predict how species will respond to potentially suitable habitats. For example, temperature and oxygen availability will interact to shape physiological and performance traits. Our model species, the wall lizard, Podarcis muralis, is a widely distributed ectotherm that continues to expand its range in Europe despite being limited by cold temperatures at high elevations and latitudes. To test the potential for organisms to expand to warming high-altitude environments, we conducted a transplant experiment to quantify the within-individual effects of high-altitude hypoxia on physiological and performance traits. Transplanted lizards maintained individual differences in physiological traits related to oxygen capacity and metabolism (hemoglobin concentration, hematocrit, and peak postexhaustion metabolic rate), as well as performance traits tied to fitness (sprint speed and running endurance). Although lizards altered blood biochemistry to increase oxygen-carrying capacity, their performance was reduced at high altitude. Furthermore, lizards at high altitude suffered a rapid loss of body condition over the 6-wk experiment, suggesting an energetic cost to hypoxia. Taken together, this demonstrates a limited potential for within-individual plasticity to facilitate colonization of novel high-altitude environments.